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Abstract 

We use food webs generated by a model to investigate the effects 
of deleting species on other species in the web and on the web as 
a whole. The model incorporates a realistic population dynamics, 
adaptive foragers and other features which allow for the construction 
of model webs which resemble empirical food webs. A large number 
of simulations were carried out to produce a substantial number of 
model webs on which deletion experiments could be performed. We 
deleted each species in four hundred distinct model webs and deter- 
mined, on average, how many species were eliminated from the web as 
a result. Typically only a small number of species became extinct; in 
no instance was the web close to collapse. Next, we examined how the 
the probability of extinction of a species depended on its relationship 
with the deleted species. This involved the exploration of the concept 
of indirect predator and prey species and the extent that the proba- 
bility of extinction depended on the trophic level of the two species. 
The effect of deletions on the web itself was studied by searching for 
keystone species, whose removal caused a major restructuring of the 
community, and also by looking at the correlation between a number 
of food web properties (number of species, linkage density, fraction of 
omnivores, degree of cycling and redundancy) and the stability of the 
web to deletions. With the exception of redundancy, we found little or 
no correlation. In particular, we found no evidence that complexity in 
terms of increased species number or links per species is destabilising. 



1 Introduction 



The world's ecosystems are increasingly being subjected to stresses that re- 
sult in large-scale changes in species population densities. These stresses 
often directly or indirectly arise from human activities and include pollu- 
tion, over-exploitation, species invasions and habitat destruction (Carlton 
and Geller 1993, Milner-Gulland and Bennett 2003). Understanding how 
ecosystems respond to such perturbations is therefore highly important. 

Here we will attempt to contribute to this understanding by focusing on 
species deletion, the complete removal of a species from an ecosystem com- 
munity, using the theoretical framework of dynamical modelling. Species 
deletion is a large-scale perturbation of particular relevance, as it is a com- 
monly used empirical tool to measure interactions strengths within real com- 
munities, and can be considered a reasonable approximation to other large 
perturbations (Paine 1980, Pimm 1980). 

The theory of small perturbations in dynamical models of ecosystems is 
well developed. It began with May's seminal work showing that the prob- 
ability of an ecosystem with random interactions being locally stable de- 
creases with both the number of species, the frequency of interactions and 
the strength of those interactions (May 1972, 1973). This result was an 
important contribution to the complexity-stability debate and contradicted 
earlier ideas that complexity should naturally lead to stability (Odum 1953, 
MacArthur 1955, Elton 1958). A crucial conceptual element to May's work 
is that only a local knowledge of the dynamics, encapsulated in the "com- 
munity matrix", is necessary to determine the stability of the population 
equilibrium to small perturbations. This is also true of perturbations that 
actually alter the position of the equilibrium, provided that they are small 
enough (Yodzis 1989). The results of small perturbations can be determined 
using only "local models" , where the species growth rates are approximated 
by linear functions of the population densities (Yodzis 2001). 

In contrast the study of large-scale perturbations requires a "global model" 
- one defined over the whole of phase space. The use of such a model will 
inevitably involve a modelling choice, but it is important that it incorporates 
phenomena such as non-linear functional responses and adaptive foraging, 
which will be likely to operate over these large changes in population density 
(Abrams 1996). 

A global model with these features has already been developed as part 
of a larger model of community coevolution by some of us (Drossel et al. 
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2001). This model assembles ecosystem communities through the repeated 
addition of new species that are modified versions of those already present. 
It is therefore conceptually similar to community assembly models (Drake 
1990, Law and Morton 1996, Morton and Law 1996, Lockwood et al. 1997), 
the difference being that new species are generated in situ rather than being 
taken from a species pool. In the model species are constantly being sub- 
jected to large perturbations in population densities as new species add and 
existing species go extinct. Crucially, species are allowed to alter their for- 
aging strategies in response to these changes. Thus the population dynamics 
of the model are particularly well suited to the study of species deletion and 
will be used in this study. Since previous studies of deletion have used the 
Lotka-Volterra equations or the equivalent discrete time Ricker dynamics, 
this will give a unique perspective on the problem (Pimm 1979, Pimm 1980, 
Borrvall et al. 2000, Lundberg et al. 2000). 

We will not only use the population dynamics of the model, we will also 
use it to generate the food webs from which species will be deleted. It 
might be preferable to use real food web structures, but then interaction 
strengths that possess a stable equilibrium for that structure would have 
to be determined, which is likely to be difficult for a large food web. The 
only studies of deletion that have used large empirical food webs circumvent 
this problem by adopting a network approach, without using an explicit 
population dynamics (Sole and Montoya 2001, Dunne et al. 2002). Their 
procedure reveals interesting features of the network topology, but by using 
the simple rule that species go extinct when none of their prey remain, they 
ignore the complex realities of the dynamics. Our model allows us to generate 
a data set of food webs with reasonably realistic topologies and interaction 
strengths (Drossel et al. 2001, Quince et al. unpublished), suitable for 
investigating the complex dynamics of deletion in multi-species communities. 

The outline of the rest of the paper is as follows. We begin with a review 
of the model in Section |2J followed by the details of the generation of the 
food web data set and the deletion experiments themselves in Section El 
The next sections concern the results of these experiments. In Section 0] we 
investigate how the trophic relationship between a pair of species influences 
the outcome of deleting one of them, in Section El we go on to consider which 
species properties correlate with large changes in community composition 
following that species removal and in Section El we consider the stability 
of the communities as a whole to deletion and the extent to which this is 
determined by their food web structure. We finish with a brief discussion of 
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the major results in Section [7| 

2 The model 

We will now give a short description of the model we will use to evolve the 
food webs to be used in our study. Further details are given in Caldarelli 
et al. (1998), Drossel et al. (2001), Quince et al. (2002) and Quince et al. 
(unpublished). 

The dynamics of the model has a different appearance depending on the 
time-scales under consideration. On time-scales of the order of the lifetimes 
of individuals, the number of species is fixed and the dynamics is determinis- 
tic: it is given by a variant of the standard equations of population dynamics. 
This includes the capacity of all species to change their predator-prey char- 
acteristics with time, competition between predators and a varying amount 
of effort that a predator puts into catching a particular prey. The system 
is allowed to change according to this dynamics for what may be relatively 
long intervals until equilibrium is reached, that is, until the populations of 
the different species present remains unchanged. If during the population 
dynamics the density of a species falls below a value N mm , usually taken to 
be 1, it is assumed to have become extinct, and is removed from the system. 

Once an equilibrium of the populations densities is reached a speciation 
event is initiated: a new species is generated by changing one of the features 
of one of the individuals of a randomly chosen species. This new species is 
then added to the system with a population density N chlld , also taken to be 1 
in these simulations. The system is then again allowed to develop under the 
deterministic equations of the population dynamics. The small population of 
the new species may give rise to a viable population, or it may die out, but 
eventually when a new equilibrium is reached, a new speciation event will take 
place. By repeating this procedure tens of thousands of times an entire food 
web can be evolved, using a combination of conventional population dynamics 
and stochastic speciation events. On these very long evolutionary time-scales 
the discrete time steps where speciation occurs are the main aspect of the 
dynamics. 

In order to define the speciation process, and more generally to be able 
to characterise a species, we need to introduce features which when taken 
together make up a species. In our model, features are specified by integers: 
a = 1, . . . , K. Any subset of L of these features constitutes a species. It is 
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assumed that the effectiveness of predator-prey relationships among species 
is due to the effectiveness of individual features against each other. Therefore 
the score of one species, i, against another, j, denoted by is defined in 
terms of the K x K matrix m a p which gives the score of feature a against 
feature [3: 



The matrix is antisymmetric. Its independent elements are random 
Gaussian variables with zero mean and unit variance chosen at the beginning 
of a simulation run and not changed during that particular run. This allows 
the score of one species against another to be calculated from ((TJ): if SV, > 
then species i is adapted for predation against species j, if SV,- = then it is 
not. We will also need to define the overlap g^, between two species i and 
j, as the fraction of features of species i that are also possessed by species 
j. The external environment is represented by a species indexed 0. This is 
assigned a random set of L features at the beginning of a run, and is not 
changed throughout the course of the run. 

Having described the structure of a species in the model, and used this to 
define the score (P) and overlap, we will now use these quantities in the con- 
struction of the population dynamics that governs the changes in population 
sizes between speciation events. The rate of change of N^t), the population 
size of species i at time t, is given by 



The function is the functional response: the rate at which one individual 
of species i consumes individuals of species j. The choice of gij essentially 
defines the nature of the population dynamics. We will give the explicit form 
chosen below, but it is here that the dependence on the score and the overlap 
functions will enter. It will also depend on the Nk(t), and so will change with 
time. 

The terms on the right-hand side of ((21) are now simply interpreted. The 
last factor represents the loss in resources for species i due to predation by 
all of the other species (in the model the measure of resources and species 
number are synonymous). The factor J2j ^i9ij on the other hand represents 
the gain to species % from predation on the set of species j including the 
environment, species 0. The environment is assigned a fixed population, 




(1) 



dNi 



N. + X^Nig^-^Njgji. 



(2) 



dt 



j j 
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Nq = Rj A, thus R is a parameter of the model that controls the rate of input 
of external resources. If it is assumed that a fraction A of the resources gained 
through predation are used to create new members of species i, the second 
term on the right-hand side of (J2J) is obtained. Finally, the first term simply 
represents the rate of death of individuals in the absence of interaction with 
other species. 

In order to briefly motivate the form of the functional response we will 
use, let us first discuss the case of a single predator i feeding on a single prey 
j. In this case 

Q (t) _ sMt) 

9tAt) ~ bN s (t) + SyNM ' [6) 

where b is a constant. We can gain more understanding of the structure of 
gij by noting that when the predators are far more numerous than the prey 
(Ni 3> Nj), Nj ~ gijNi: the feeding rate of the predators is limited only 
by the number of prey. In the other limit, when the prey is very abundant 
compared with the predators (Nj 3> N), ~ Sij/b: each predator feeds 
at a constant maximum rate. This latter result also gives an interpretation 
to the constant b. Having introduced the basic form (jHJ), we can now state 
general form for the functional response used in the model: 

(f\ Sjjfij(t)Nj(t) 

9ij{ ] bN 3 {t) + E fe a ki S kj f kj (t)N k (t) • 1 ] 

There are two new aspects present in (@J) and absent in Q: 

1. Interference competition between predators of prey species j is mod- 
elled by the factor a ki . We take a„ = 1 and a ki < 1, i k to reflect the 
fact that competition between members of the same species is typically 
stronger than competition between different species. In fact, we expect 
that the more species are alike, the greater will be the competition 
between them, and therefore take 

ay = c + (1 - c)qij , (5) 

where c is a constant lying between zero and one which is the residual 
degree of competition that exists even if two competing predators have 
no features in common. 

2. Adaptive foraging is modelled using the factors /y. The effort can 
be viewed as the fraction of time an individual of species i spends 
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predating on species j or the fraction of the population of species i 
dedicated to consuming only j. These efforts must satisfy J2j fij — 1 
for all i. To determine the it seems reasonable to assume that the 
gain which an individual of species i makes in consuming individuals 
of species j (that is, gr^), divided by the amount of effort i puts into 
this task (that is, fy), should be the same for all prey species j. Using 
this condition, together with the normalisation of the /y, leads to 

fij(t) = J ij{t ] r ■ (6) 
J2k 9ik(t) 

This choice of efforts can be shown to be an evolutionary stable strategy 
(Drossel et al. 2000), or in the terminology of foraging theory, an ideal 
free distribution of predators across prey (Fretwell and Lucas 1970). 

The calculation of the efforts, through (JHJ), effectively introduces a new 
behavioural timescale into the problem. We assume that the efforts change 
on a much shorter timescale than the population densities Nj(t), and there- 
fore that they may be found by iterating (jlj) and (jOJ assuming constant 
population densities. When this process has been completed, we may then 
move on to updating the population densities. Thus we do not treat the 
efforts as dynamical variables, instead we assume that they are a function of 
the population densities, even if we have no explicit form for this function. 
This describes the essential features of the model. In the rest of the paper 
we will explore the consequences of performing deletion experiments on food 
webs created by the model. 



3 Deletion experiments 

The set of model food webs used in this study was obtained by performing 
four hundred independent simulations of the model, each simulation lasting 
for 120000 speciations. These simulations were independent, in that different 
pseudo-random number sequences were used in their generation. Thus the 
simulations differed in their random matrices, m a p, environment features and 
speciation events. The same parameters were used in all the simulations these 
being R=lx 10 5 , b = 5 x 10~ 2 , c = 0.5, A = 0.1 and N m[n = N child = 1.0. 
The effect of altering the model parameters is investigated in Quince et al. 
(unpublished). The four hundred final food webs from these simulations, 
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Figure 1: A typical model food web. 



which after 120000 speciations will have structures drawn from a stationary 
distribution, constituted the ecosystem data set. 

A typical food web is shown in Fig. Q Each species in this diagram is 
represented as a circle, the sizes of which are the same for all species. This 
differs to the convention adopted in Quince et al. (2002), where the radii 
of the circles were proportional to the logarithm of the population densities. 
The arrows represent predator-prey interactions, with the arrow pointing 
from the prey to the predator. The intensity of the arrow is proportional to 
the fraction of the predators diet that consists of that particular prey. The 
vertical arrows originating from the base of the diagram, rather than from 
another species, indicate that the species is feeding off the environment. 
The species are positioned vertically according to trophic height, defined as 
the average path length from the species to environment, the average being 
weighted by predator diet fractions. The dashed lines show the position of 
integer values of the trophic height. 

The four hundred food webs in the model data set span a range of struc- 
tures. We shall quantify this variation with five food web statistics: 

1. The total number of species in the food web denoted by S. 

2. The links per species denoted by L/S. This quantity is simply the 
number of predator-prey interactions divided by the number of species, 
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Statistic 


Symbol 


Mean 


Std. Dev. 


Number of species 


S 


63.8275 


7.5152 


Links per species 


L/S 


1.6881 


0.1258 


Fraction of omnivores 





0.1606 


0.0474 


Degree of cycling 


c 


0.0049 


0.0022 


Redundancy 


V 


0.1884 


0.0706 



Table 1: The means and standard deviations of five food web statistics for 
the four hundred model food webs in the data set. 

where we will use the convention of counting a link if it constitutes 
greater than 1% of a predator's diet. 

3. The fraction of omnivorous species in the food web denoted by O. We 
define omnivorous species as those which feed at more than one trophic 
level and define the trophic level of a species to be the shortest path 
from that species to the environment. The justification for this choice is 
the observation (Yodzis 1984) that the shortest path between a species 
and the environment tends to be the most important energetically. 

4. The degree of cycling denoted by C. The method we adopt for mea- 
suring the amount of cycling of energy in the food webs is based on 
both the ideas presented in Ulanowicz (1983) and on source code kindly 
provided by the author. Essentially a backtracking algorithm was first 
used to identify each cycle in a food web. Having done this, the amount 
of energy flowing in a cycle was identified with the strength of the weak- 
est link of the cycle, exactly as in Ulanowicz (1983). To define C we 
measure the proportion of a predator's prey obtained through cyclic 
flow averaged over the species in the web. 

5. The ecosystem redundancy denoted by V. The redundancy of an 
ecosystem is the proportion of species which can be considered super- 
fluous to the functioning of the ecosystem as a whole (Walker 1992). 
We define a species to be redundant in our food webs if at least one 
other species possesses the same pattern of trophic links i.e. the same 
predators and prey. As for the calculation of L/S, only links form- 
ing greater than 1% of the predator's diet are used in this calculation. 
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Then the redundancy, V, is the fraction of redundant species in the 
food web. 

In Table ^the means and standard deviations of the five statistics are shown 
for the four hundred food webs in the data set. 

This data set was then used to investigate the effect of species deletion. A 
single species was removed from a web and the population dynamics iterated 
until a stable equilibrium was reached. If, during this process, the population 
of a species fell below iV mm , then it was removed ("went extinct") in accor- 
dance with the criterion applied when evolving the communities. For each 
web every species was deleted independently, that is, the webs were returned 
to their original state between deletions. In total, over all the food webs, 
25531 species were deleted. In Fig. |2]we show the frequency distribution of 
number of further extinctions for these deletions. This distribution does not 
decay exactly exponentially with further extinction number, but it does have 
a characteristic size of just a few species. The number of further extinctions 
seems bounded, the largest number is seventeen, and on no occasion is the 
whole web of typically sixty species close to collapse. The effect of deleting 
a species is localised in the webs. 

4 Trophic relationships and species deletion 

In this section we investigate how the trophic relationship between a pair 
of species influences the impact that deleting one of the pair will have on 
the other. We began by simply taking the four hundred webs described in 
the previous section and placing each species into four non-exclusive cate- 
gories according to its trophic relationship with the deleted species. These 
categories were predators, prey, competitors and indirect predators. The def- 
initions of the first two categories are obvious. A "competitor" was defined 
as any species that shared a prey with the deleted species. The final category 
of "indirect predators" requires a little more explanation. Consider any two 
species in an ecosystem, i and j. If there exists one or more paths from j to 
i travelling only up trophic links, from prey to predator, and the minimum 
length of those paths is greater than one, then we define species i to be an 
indirect predator of species j. In other words, species % is an indirect preda- 
tor of j if some of the resources consumed by i came originally from j, but 
species i does not prey on j directly. In calculating the categories all links 
which formed less than one percent of the predators diet were ignored. 
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Figure 2: The frequency distribution of the number of further extinctions for 
the 25531 species deletions. 

Having defined these categories we calculated two quantities. The first 
was the mean number of species in each category that went extinct following 
species deletion, with the average performed over every deleted species in 
every web. The second was the probability that a species with a particular 
trophic relationship to the deleted species went extinct. This quantity was 
estimated by simply averaging the proportion of species with a particular 
relationship that went extinct. These results are shown in Table El The first 
row shows the results for all species, in order to aid comparison with the other 
results. We see that, on average, following the deletion of a species from a 
web, a further 1.366 species went extinct and that any given species had a 
probability of 0.021 of going extinct. A clearer picture is usually obtained 
by using probabilities, rather than mean values, when making comparisons 
between tropic relationships. This is due to the weaker dependence that 
probabilities have on the number of species with that particular relationship 
to the deleted species. From these results we see that that a predator of a 
deleted species is more than ten times as likely to go extinct following the 
deletion of its prey than an average species in the web. This makes ob- 
vious intuitive sense. More interestingly it is also found that a prey of the 
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Trophic relationship 


Mean number of extinctions 


Probability of extinction 


All species 


1.366 


0.021 


Predators 


0.675 


0.233 


Prey 


0.132 


0.041 


Competitors 


0.013 


0.004 


Indirect predators 


0.658 


0.106 



Table 2: The effect of trophic relationship on the probability and mean num- 
ber of species going extinct following deletion of a species from an ecosystem. 

deleted species is twice as likely to go extinct as an average species in the web 
when its predator is removed. The most likely explanation for this is that 
an effect known as predator-mediated coexistence or keystone predation is 
operating. This is a mechanism whereby a predator allows inferior competi- 
tors to coexist with a superior competitor by predation of the competitively 
dominant species. Thus the deletion of the predator can lead to extinctions 
amongst its prey. This has been observed in real communities (Paine 1974, 
Lubchenco 1978, Navarrete and Menge 1996) and in theoretical studies of 
simple dynamical systems (Fujii 1977, Shigesada and Kawasaki 1988). 

The plausibility of the above mechanism is supported by the fact that 
competitors of the deleted species are five times less likely to go extinct 
than an average species, suggesting that competition plays an important 
part in structuring these communities. The effect on indirect predators is also 
significant: they are five times as likely to go extinct following the removal 
of their indirect prey compared to the average. It is interesting to note that 
the effect on an indirect predator is less than that on a direct predator. We 
will return to this later. 

4.1 Predator-prey interaction strengths and extinction 
probabilities 

The results in Table |2] show that removing the prey, or indirect prey, of a 
species has a negative effect on the predator. Clearly the importance of a 
prey species to its predator will vary between predator-prey pairs. A measure 
of the positive effect that a prey species has on the predator will be the effort 
fij, corresponding to the fraction of the prey in the predators diet (Ulanowicz 
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and Puccia 1990). 

We can derive a similar quantity for the positive effect that an indirect 
prey species has on its indirect predator. Consider the square of the / matrix 
fij — J2k=o fikfkj, where the sum is over all species in the web and k — 
corresponds to the environment. Since fkj is the fraction of species fc's diet 
that comes from species j, then fikfkj is the fraction of species i J s diet that 
comes from species j via species k. If we now sum k over all species in 
the web, we obtain the fraction of species z's diet that comes from species j 
through all paths of length 2. Therefore if we define the matrix F = Y^=i / n > 
then is the fraction of species i's diet that comes from species j via all 
possible paths. Since all resources originally derive from the environment, 
F iQ = 1 for any species i. For the direct matrix f\j, normalisation ensures 
that Ylj=ofij = 1-0- However for the indirect matrix, 



with equality only when species i is a basal species. 

If the infinite series of matrices converges, it can be calculated using 
F — (1 — f) f. All eighty / matrices in this study were such that the sum 
was finite, however we have not been able to prove the convergence of the 
sum for a general web. 

We are interested in separating the effects of indirect and direct predation. 
This can be achieved by defining a further matrix I such that 



Thus if Iij > 0, then species % will be an indirect predator of species j 
according to the definition given above. 

We can now examine how the probability that a predator or indirect 
predator goes extinct, following removal of its prey, varies with the proportion 
of the prey in the predators diet. This was carried out by dividing all the 
predator-prey (indirect predator-prey) pairs in the four hundred ecosystems 
into bins of size 0.05 according to their associated fij (iy) values. The 
proportion of prey removals that resulted in the extinction of the predator 
(indirect predator) was then calculated for each bin and plotted in Fig. El 

The interpretation of Fig. |3] is quite straightforward. Ignoring fluctua- 
tions attributable to sampling effects, both curves are monotonic. As the 
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Proportion of prey in predator's diet (f or I..) 

Figure 3: The fraction of predators (indirect predators) going extinct upon 
prey removal as a function of the proportion of that prey (indirect prey) in 
the predators diet. 

fraction of the predators, or indirect predators, diet that is obtained from 
the deleted species increases, so does the probability of predator extinction. 
However the probability of extinction does not become large until the prey 
constitutes a significant fraction of the predators diet. For direct predators 
it does not reach 10% until ~ 0.65 and for indirect predators this occurs 
when Iij « 0.45. This probably arises from incorporating adaptive foraging 
into the population dynamics: predators can survive events that remove a 
large portion of their prey. In both cases the extinction probability rapidly 
approaches, but does not quite reach, 100% as (Jy) approaches 1.0. In 
fact, when a prey constitutes greater than 99% of a predators diet its removal 
leads to predator extinction only 97% of the time. This again illustrates the 
effect of adaptive foraging: 3% of the time the predator must be surviving 
by exploiting a prey species that previously formed less than 1% of its diet. 

The two curves in Fig. |3] have a similar form, but at any given or 
Iij value the probability of extinction is greater for indirect predators. This 
seems to contradict Table 121 which shows a greater probability of deletion for 
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Figure 4: The probability distribution of and for all > 0.01 and 
Jjj > 0.01. These results are a compilation over the four hundred food webs 
described in Sectional 

direct predators than indirect predators. However this can be explained by 
the distribution of the non-zero elements of the I and / matrices. If there are 
more small values of Jjj than fy, then defining a predator (indirect predator) 
as having > 0.01 (Jjj > 0.01) will lead to a lower probability of deletion 
for indirect predators than direct predators, even though for a given value of 
fij (A?) the e ff ec t °f removing the indirect predator is more significant. That 
this is the case can be seen from Fig. EJ where the probability distributions 
of and Iij for all > 0.01 and 1^ > 0.01 are shown. This does not of 
course explain why an indirect predator is more likely to go extinct when 
prey constituting a given fraction of its diet is removed. This may arise from 
the assumption that the net importance of multiple paths to the predator can 
be obtained by simply summing their individual weights implicit in Eq. (J7J). 

We showed in Table |2] that a prey species has an increased probability 
of extinction following removal of its predator. We might also expect that 
this effect will depend on the fraction of the prey in the predators diet. This 
is investigated in Fig. where the probability of prey extinction following 
predator removal has been estimated as a function of f\j, by placing all the 
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Figure 5: The fraction of prey going extinct upon predator removal as a 
function of the proportion of that prey in the predators diet. 

predator-prey interactions from the four hundred food webs into bins of size 
0.1, and calculating the proportion of deletions that led to extinction of the 
prey. From this graph we see that prey extinction probability peaks at an 
intermediate value around = 0.25, where slightly less than 15% of prey 
are going extinct. This is what we would expect if the mechanism for prey 
extinction is predator-mediated coexistence as proposed above, since if a prey 
forms only a fraction of a predators diet, then that predator is also likely to 
be exploiting its competitors. 

4.2 Trophic Levels 

In TableEJwe categorised species according to their trophic relationship to the 
deleted species. Another way to categorise species in a food web is by trophic 
level. As mentioned in Section 01 we use the minimum path length definition 
of trophic level (Yodzis 1984). Having assigned a trophic level to each species 
in our ecosystems, we calculated the mean number of species on each level 
that went extinct as a function of the trophic level of the deleted species. 
These results are shown in Table HU where the same averaging procedure was 
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Level 


All 


1 


2 


3 


4 


All 


1.366 


0.011 


0.649 


0.679 


0.026 


1 


6.571 


0.012 


3.651 


2.783 


0.124 


2 


0.705 


0.013 


0.121 


0.560 


0.011 


3 


0.544 


0.009 


0.359 


0.161 


0.015 


4 


1.223 


0.006 


0.498 


0.709 


0.010 



Table 3: The mean number of species on a given level that go extinct as a 
function of the trophic level of the deleted species. The rows refer to the 
trophic level of the deleted species and the columns to the trophic level in 
which further extinctions occurred. 

used as for Table that is, the average was performed over all the deleted 
species in all eighty ecosystems. 

The analysis of this data is aided by statistics on how the number of 
species S, the average prey number S prey , and average predator number 
Spredator, var y between trophic levels. The number of species S in each trophic 
level was calculated by averaging over the eighty webs in this study. The other 
statistics were averaged over the total number of species in each trophic level. 
These statistics together with standard deviations in brackets are shown in 
Table H 

Several patterns emerge from TableEJ For example, if we remove a species 
from a particular level then this leads to extinctions on the level above and 
this effect diminishes as we increase the level of the deleted species. Thus 
the mean number of level 2 species going extinct upon removal of a level 1 
species is 3.651, which is roughly seven times the mean number of level 3 
species going extinct following a deletion on level 2 (0.560), which in turn is 
about forty times larger than the mean number of extinctions in level 4 as 
a result of a deletion on level 3 (0.015). We can explain this with reference 
to Table EJ from which we see that while the mean number of prey increases 
rapidly with trophic level, the mean number of predators decreases. These 
two effects will be complementary, so that as we increase the trophic level 
of the deleted species, the number of predators of that species decreases, 
and those predators become less specialised on the deleted species. This 
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Level 


S 


<7 

"-'prey 


'S'predator 


1 


7.76(1.06) 


0(0) 


4.26(1.34) 


2 


30.19(4.41) 


1.13(0.37) 


1.68(0.74) 


3 


24.63(4.27) 


2.71(1.62) 


1.25(0.49) 


4 


1.26(0.84) 


5.84(3.31) 


1.09(0.29) 



Table 4: The dependence of number of species S, number of prey Sp Te yj and 
number of predators S^dator, on trophic level. The figures in brackets are 
standard deviations. 

leads to less extinctions in the level above when the species is removed. As a 
concrete example, consider deletion of a level 1 species. This will, on average, 
affect about 4 species on level 2 and these 4 species will be almost entirely 
dependent on the deleted species (as mentioned previously, deletion of the 
only prey of predator leads to extinction of that predator almost 100% of 
time). Thus we would expect about 4 predators in level 2 to go extinct as a 
result of a deletion on level 1. This compares well with the actual value of 
3.651. We also see from Table El that the effect of deleting a species on the 
level above propagates up trophic levels, so that removing a level 1 species 
leads to deletions on level 2, which in turn leads to extinction on levels 3 and 
4. This is as we would expect, given that we have already illustrated that 
consequences of species deletion can propagate up food webs. 

As well as extinctions occurring in the levels above a deleted species, 
we find that deleting a species can lead to extinctions in the level below 
it. The magnitude of this effect seems to decrease with the trophic level of 
the deleted species: the mean number of level 3 species going extinct given 
deletion on level 4 is 0.709, compared to 0.359 for the number of level 2 
species going extinct as the result of a deletion in level 3. Similarly, the 
deletion of a level 2 species leads to an average of just 0.013 extinctions on 
level 1. These observations can be explained if we assume that predator- 
mediated coexistence causes the extinctions, since the more prey a predator 
has, the stronger we might expect this effect to be. Then the pattern of 
increasing extinctions in the level below, with increasing level of the deleted 
species, merely reflects the increase in average prey number with trophic 
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level observed in Fig. The question of why the model food webs show 
these changes in mean predator and prey number with trophic level is not 
addressed here, but is discussed in Quince et al. (unpublished). 

5 Keystone species 

Thus far the focus of this study has been on pairs of species and how the 
trophic relationship between them influences the effect that deleting one 
species will have on the other. It was found that the effects of deletion 
can propagate both up and down food webs. We now consider the related 
question of what factors determine the impact that deleting a species will 
have on the whole web. Specifically, we ask whether there are any consistent 
differences between species whose removal causes little change in the food 
web, and those which play a major role in structuring the community. The 
latter are sometimes referred to as "keystone species", although this term 
can be restricted to those species whose importance is large relative to their 
population size (Power et al. 1996). Here we will simply be interested in iden- 
tifying factors that statistically influence the number of further extinctions 
that follow the removal of a species. 

The trophic level of the deleted species has an effect on the expected 
number of further extinctions. This is shown in the first column of Table E3 
where the mean number of further extinctions in the web as a whole is seen 
to be largest for species on trophic level 1. This probably reflects the greater 
importance of effects propagating up, rather than down, the food web. An 
alternative way to categorise species according to trophic role, which is better 
at separating these two processes, is into the following three classes: 'top' 
(species with no predators), 'intermediate' (species with both predators and 
prey) and 'basal' (species with no prey). When this is done, it is found that 
removal of a basal species causes an average of 6.58 further extinctions as 
opposed 0.65 and 0.61 for intermediate and top species respectively. Com- 
paring the figures for basal and top species reveals that top-down processes 
are indeed statistically less important than bottom-up in these communities. 

The above categories are quite broad. One property which we can use to 
more finely discriminate between species is the number of other species they 
interact with through predator-prey links. We might expect this to correlate 
with the impact of deleting the species on the food web. That this is indeed 
the case is shown in Fig. |Hl where the average number of further extinctions 
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Figure 6: The number of further extinctions following the deletion of a species 
as a function of its total number of predators and prey also known as the 
node degree k. The results are subdivided according to whether the deleted 
species is basal, intermediate or top and are averaged over all species in all 
four hundred communities. The error bars give standard errors in the mean. 

is plotted as a function of the node degree k in network terminology. In the 
case of food webs the latter corresponds to the number of predators plus the 
number of a prey of a species. The results are shown for basal, intermediate 
and top species separately. For all three classes the average number of further 
extinctions increases with k, so the most connected species are the ones whose 
removal has the greatest effect on the food web structure. The results for 
top species provide further support, albeit circumstantial, that the top-down 
effect in our webs is predator-mediated coexistence, since we would expect 
the importance of this effect to increase with the number of prey of the 
predator. They can be contrasted with an earlier study that failed to find 
such a relationship (Pimm 1980). 

These results have relevance for the studies of real food web robustness 
to deletion mentioned in the introduction. These obviated the need for an 
explicit dynamics by considering multiple removals, and judging a further 
extinction to have occurred when all the prey of a species are absent (Sole 
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Figure 7: The number of further extinctions following the deletion of a basal 
species % plotted against its bottom- up keystone index (K l b ). The faint grey 
dots give individual data points present in order to give a sense of the dis- 
tribution. The black diamonds give means with error bars showing standard 
deviations of the distribution of K\ values for each number of further extinc- 
tions. 

and Montoya 2001, Dunne et al. 2002). These studies found that removing 
the most connected species resulted in more secondary extinctions for the 
same number of species removed. Our results suggest that if a population 
dynamics was included in these studies, the food web structures would be 
even more sensitive to the removal of highly connected species. 

The correlation between k and further extinction number is very strong, 
especially given that it quantifies the number of direct interactions of a 
species, and does not give information on the strength of indirect effects, 
which we have already shown to be important (Jordan and Scheuring 2002). 
It is difficult to devise measures of species importance that do incorporate 
indirect effects. Potential candidates are the bottom-up and top-down key- 
stone indices of Jordan et al. (1999). These were originally devised for binary 
food webs, but is easy to extend them to the model food webs including diet 
compositions considered here. In fact the bottom-up keystone index for a 
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Figure 8: The number of further extinctions following the deletion of a top 
species % plotted against its top-down keystone index (Kl). The faint grey 
dots give individual data points and the black diamonds give means with 
error bars showing standard deviations of the distribution of K\ values for 
each number of further extinctions. 

species i is simply 

K l = J2 F ii- ( 9 ) 

3=1 

Here the matrix F is defined by F = Y£?=i /" an d, as discussed in the 
previous section, its elements give the fraction of species z's diet that comes 
from species j via all possible paths. This makes the meaning of Kl clear: 
it is the total number of species that depend on i directly or indirectly for 
resources. The complementary quantity K\ measures the strength of top- 
down effects it can be defined in a similar way for non-binary food webs, which 
we will not discuss here. There is in fact a very good correlation between K l h 
and further extinction number for basal species. This shown in Fig. [7| This 
suggests that this quantity is effective at predicting the strength of bottom- 
up effects. We also found a correlation between K\ and the number of further 
extinctions following the removal of a top species (Fig. |HJ) although not as 
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strong as for basal species. Furthermore the composite quantity K l = K\+K l b 
did have some success at predicting the outcome of removing an intermediate 
species. This suggests that although top-down effects might be more difficult 
to quantify with these network measures, and they are not yet as effective 
as simply counting the number of predators and prey, they may well prove a 
practical way of predicting the importance of species in food webs. 

6 Stability of food webs to deletion 

In the previous section we focused on the effect that deleting one species has 
on another in the model food webs. Here we will consider the stability of 
the communities as a whole. In particular we will ask whether that stability 
is determined by the topological structure of the ecosystem. We begin by 
defining deletion stability as the fraction of species that have the property 
that, when they alone are deleted, further extinctions do not take place. This 
definition is similar to that given by Pimm (1979), and has the advantage that 
it should not have any in-built dependence on species number. Its calculation 
is illustrated in Fig. 

We will use the four hundred model food webs detailed in Section |3] as 
our data set and describe their topological structure with the five food web 
statistics also defined there. These five statistics quantify four food web prop- 
erties that it has been suggested may impact ecosystem stability: complexity 
in terms of number of species and linkage density (S, L/S), the amount of 
omnivory in the web (O), the importance of cycles (C) and the redundancy 
(V). We remind the reader that the means and standard deviations of these 
statistics for the data set are given in Table ^ 

If we possessed a range of food webs that vary for each property indepen- 
dently, whilst the other properties remain fixed, then this analysis would be 
quite simple. Indeed this was the approach adopted in previous studies us- 
ing small pre-defined food web structures (Pimm 1979, Pimm 1980, Borrvall 
et al. 2000). However by generating the food webs by evolving a number 
of communities at a particular set of parameter values, we are faced with a 
less straightforward situation. The food webs vary for all properties simul- 
taneously, and some of the properties are significantly correlated with one 
another. In fact, the situation is much like analysing real food web struc- 
tures. The advantage is that our range of structures are much more complex 
and realistic, the disadvantage is that we are forced to adopt a statistical 
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Figure 9: The food web of Fig. [T] with those species whose deletion causes 
further extinctions shaded grey. There are 27 such species. The number of 
species in the web is 61, giving a deletion stability, Sd = 34/61 or 0.557. The 
latter is a fairly typical value for the webs in this study. 

approach. 

We will commence our analysis by examining the correlation matrix for 
all the variables, both the dependent variable Sj, and food web properties. 
This is shown in Table |5] The deletion stability is significantly correlated 
with the fraction of omnivores (negatively) and the redundancy (positively). 
However because of the many significant correlations between the variables 
themselves, we can not conclude that a smaller proportion of omnivores or 
greater redundancy will be associated with higher deletion stability, all other 
properties being unchanged. 

One way to shed some light on this problem is to use multivariate linear 
regression (Jobson 1991). This statistical procedure essentially assumes that 
the dependent variable, in this case Sd, can be described as a linear function of 
the independent variables, the five food web properties plus a noise term. The 
results of such an analysis are shown in Table El Examining the regression 
coefficients, bj, for the individual variables we see that two, the redundancy 
again and the links per species, are deemed to have a strong influence on the 
deletion stability, and that both of these effects are positive. More pertinently 
however, examining the fit as whole we see from the 1Z 2 value that only 8% 
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Sd 


s 


L/S 





C 


V 


s d 




-0.054 


0.063 


-0.162* 


-0.111 


0.197* 


s 






0.354* 


0.073 


-0.129 


0.034 


L/S 








0.197* 


-0.196* 


-0.264* 













0.359* 


-0.282* 


c 












-0.085 



Table 5: The correlation matrix for the variables: deletion stability (Sd), 
number of species (S), links per species (L/S), fraction of omnivores (O), 
degree of cycling (C) and redundancy (V). The values given are the linear 
correlation coefficients r between the pairs of variables calculated over the 400 
food webs. The correlations judged to be significant, those with a probability 
p of no correlation smaller than the Bonferroni corrected value of 0.05/15, 
are highlighted with an asterisk. 

of the variation in stability is being explained by these variables. 

We therefore conclude that for this data set a fairly robust positive rela- 
tionship exists between the proportion of redundant species in the web and 
the stability of the web to deletion. This result makes sense. We defined 
a redundant species as one that is functionally equivalent, in the sense of 
possessing the same predators and prey, to at least one other species. Thus 
the removal of a redundant species is unlikely to cause further species to 
go extinct, since its functional equivalents should be able to increase their 
population sizes and compensate for the loss. This supports the hypothesis 
that increased redundancy in ecosystems will result in increased functional 
reliability (Walker 1992, Naeem 1998). 

We find no evidence that complexity in terms of increased species number 
or links per species is destabilising. In fact there was evidence that increasing 
the latter actually reduced the probability of further extinctions when species 
are deleted. This can be compared to early work on deletion stability in 
small food webs modelled with Lotka-Volterra dynamics, where it was found 
that increasing either the number of species or connectance, L/S 2 , rapidly 
decreased stability (Pimm 1980). The difference can probably be attributed 
to the more realistic structures and global dynamics used here, in particular 
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Variable 


bj 


Std. err. 


t 


V 


Intercept** 


0.4157 


0.0803 


5.1766 


0.0000 


S* 


-0.0016 


0.0007 


-2.2903 


0.0225 


L/S** 


0.1436 


0.0453 


3.1738 


0.0016 


0* 


-0.2512 


0.1183 


-2.1231 


0.0344 


c 


-1.4098 


2.4820 


-0.5680 


0.5703 




0.3043 


0.0749 


4.0634 


0.0001 



Table 6: The multivariate linear regression of deletion stability (Sd) as a func- 
tion of number of species (S), links per species (L/S), fraction of omnivores 
(O), degree of cycling (C) and redundancy (V). This fit had 1Z 2 = 0.08186, 
T = 7.026 on 5 and 394 degrees of freedom, and P = 2.636 x 10~ 06 . The co- 
efficients that are significantly different from zero (p < 0.05) are marked with 
an asterisk, and those that are highly significant (p < 0.01) with a double 
asterisk. The calculation was performed using the software package S-Plus 
6.0. (Mathsoft Inc. 2000) 

incorporating adaptive foraging into the population dynamics. 

This result can also be compared to more recent work examining the 
effect of deleting species from simple three level food webs, constructed such 
that all species on a trophic level were functionally equivalent (Borrvall et 
al. 2000). These authors found that stability to deletion increased with the 
number of species on each trophic level. Their set-up corresponds to keeping 
V constant and equal to one, whilst increasing the total number of species, S, 
a variable for which no effect was observed here. This study also used Lotka- 
Volterra dynamics and this may explain the discrepancy. In particular, they 
observed that removing a predator had no effect on its prey, in contrast to 
the results detailed in Section E| 

We also found that omnivory and cycling were unimportant in determin- 
ing deletion stability. In fact, all our structural properties taken together 
explained very little of the variation between the webs. It may be that the 
measure of deletion stability we used is intrinsically noisy — it does seem 
sensitive to the presence of one or two vulnerable species — or it may be 
that our statistics are not capturing the properties that are important in 
determining robustness to deletion. 
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7 Conclusion 



In this paper we have shown that deletion experiments, which are very dif- 
ficult and time-consuming to carry out in real communities, can be easily 
implemented on model webs, which previous studies have shown have many 
of the characteristics of real webs (Drossel et al. 2001). We found a num- 
ber of interesting results on species removal from communities which differed 
from some previous studies. The food webs as a whole were shown to be 
quite robust to deletions, and individual species were able to survive the loss 
of prey constituting a major fraction of their diet. Deletions were shown to 
cause further extinctions amongst species both above and below the deleted 
species in the food web. These phenomena arose out of the complex popula- 
tion dynamics used in the model which incorporates adaptive foraging. They 
illustrate the importance of using a realistic global dynamics when consider- 
ing community responses to large scale perturbations such as deletion. This 
contrasts with studies that either lack an explicit dynamics (Sole and Mon- 
toya 2001, Dunne et al. 2000) or use Lotka-Volterra equations (Pimm 1979, 
Pimm 1980, Borrvall et al. 2000, Lundberg et al. 2000). 

In addition to a realistic dynamics, effective studies of species removal 
require a realistic set of structures. By using a set of large webs composed of 
species spanning a range of trophic roles, we were able to show that removing 
the most connected species resulted in the most further extinctions and that 
recently developed 'keystone indices' were fairly effective at predicting species 
importance. A study using small simple webs failed to find these relationships 
(Jordan et al. 2002). The range of food web structures studied allowed us 
to show the important role that redundant species play in increasing food 
web robustness to deletion, and find that there was no correlation between 
increased complexity and decreased stability to deletion. This adds another 
component to the stability-complexity debate. 

This paper represents the first attempt to join complex food web struc- 
tures with realistic population dynamics to study species loss from communi- 
ties. It would be interesting to see if the results would change if real, rather 
than evolved, food web structures were used with such a dynamics or if the 
particular choice of dynamics were changed. In any case, we believe that 
we have shown that it is crucial for models to display a degree of realism, 
if reliable deductions concerning the consequences of species deletions are to 
be made. 
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